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MEMORANDUM  REPORT  NO.  953 


JMSchmidt/ sdb 

Aberdeen  Proving  Ground,  Maryland 
December  1955 


ANALOG  COMPUTER  SIMULATION  OF  FLIGHT  CHARACTERISTICS 
OF  TWO  90MM  FIN  STABILIZED  SHELL  T108E40  AND  T316 

ABSTRACT 

A  study  was  made,  with  the  aid  of  an  analog  computer,  of  the  effect 
of  various  spin  histories  upon  the  dispersion  of  two  90mm  shell,  the  T108 
and  T316.  The  effects  of  resonance  are  discussed  and  graphs  illustrating 
the  motions  resulting  from  several  different  rates  of  spin  are  presented. 
A  comparison  is  made  between  the  flight  characteristics  of  the  two  shell. 
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damping  rates  of  epicyclic  arms 
turning  rates  of  epicyclic  arms 
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INTRODUCTION 


Recent  firings  of  sample  groups  of  the  T108E40  shell  (Figure  l)  in 
current  production  indicated  extremely  poor  accuracy  in  comparison  with 
the  accuracy  of  the  same  shell  during  its  developmental  stages.1  There 
naturally  aiose  the  problem  of  determining  the  causes  of  this  incongruity . 
The  current  production  shell  does  not  differ  significantly  in  its  physical 
characteristics  from  its  prototype.  However,  the  tail  which,  during 
development,  was  made  from  aluminum  73ST,  is  now  made  from  aluminum  63ST, 
which  is  known  to  be  weaker.  Experimental  firings  in  the  Transonic  Range 
showed  that  the  63ST  tail  very  frequently  deforms  during  launching.  These 
firings  also  revealed,  moreover,  that  as  a  result  of  tail  deformation 
the  axial  spin  and  subsequent  spin  history  of  this  shell  in  flight  varies 
over  a  very  wide  range. 

Unfortunately,  the  nature  of  this  shell  is  such  that  it  must  remain 
within  rather  narrow  bounds  of  spin.  The  upper  bound  is  imposed  by  penetra¬ 
tion  requirements  and  is  about  25  RPS ;  the  lower  bound  is  the  zone  of 
resonance  at  about  8-2  RPS.  The  initial  spin  is  imparted  to  the  shell 
by  the  plastic  obturating  ring  and,  depending  upon  the  characteristics  of 
this  ring,  usually  lies  between  12  and  15  RPS.  Normally,  the  initial  spin 
of  1.9°/ft.  will  decay  to  about  l°/ft.  at  2000  yards,  a  condition  which 
could  be  described  as  satisfactory.  In  the  case  of  a  shell  with  a  deformed 
tail  the  normal  spin  decay  might  well  be  accelerated  and  the  shell  would 
approach  the  resonance  condition  early  in  its  flight,  the  resonant  spin 
for  this  shell  being  l.l2°/ft.  It  was  felt  that  the  build-up  in  the  size 
of  the  yawing  motion  which  characterizes  the  resonance  phenomenon  could 
have  a  detrimental  effect  upon  the  firing  accuracy  of  the  missile. 

In  connection  with  this  problem  a  new  development  came  into  promi¬ 
nence,  namely  the  9Ctam  T3l6  shell  (Figure  2).  The  T3l6  is  a  shorter  shell 
with  a  spiked  nose  with  the  result  that,  at  the  expense  of  increased  drag, 
the  destabilizing  lift  on  the  body  of  the  shell  is  greatly  reduced  and  the 
center  of  pressure  of  the  body  is  moved  closer  to  the  center  of  gravity 
of  the  missile.  Therefore,  to  be  stable,  this  shell  requires  a  relatively 
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short  boom.  It  has  no  obturating  ring  and  its  initial  spin  is  much  lower 
than  that  of  the  T108.  Since  the  critical  frequency  of  the  TJl6  is 
2.54°/ft. ,  its  spin  must  increase  from  its  initial  zero  or  very  low  value 
to  a  relatively  high  rate  of  spin  in  order  to  get  into  resonance.  The 
likelihood  of  this  occurrence  is  diminished  by  the  fact  that,  since  the 
shell  is  quite  compact,  there  is  less  opportunity  for  the  tail  to  be 
damaged  as  the  shell  leaves  the  gun. 

This  report  deals  with  the  investigation,  by  means  of  analog  computer 
simulation,  of  the  effects  of  various  spin  histories  on  the  trajectories 
of  these  two  shell.  Since  resonance  appeared  to  be  a  possible  reason  for 
the  poor  firing  accuracy  of  the  T108E40  shell,  particular  attention  is 
given  to  resonance  conditions.  To  this  end,  six  spin  histories  were 
arbitrarily  assigned  to  each  shell,  as  shown  in  the  following  table. 
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TABLE  1 


SPIN  HISTORIES  STUDIED 

* 

Resonance 


Case 

(d©/dz)Q 

(dO/dz)^ 

Occurs  at 

(deg. /ft. ) 

(deg. /ft, ) 

(yards) 

T108E40 

IA 

0 

0 

1,2 

1.9 

1.9 

1,3 

1.9 

0 

(1601)** 

1.9 

(-3.2) 

500 

1,5 

1.9 

(-0,9) 

1000 

1,6 

1.9 

(-0.01) 

1500 

T316 

11,1 

0 

0 

— 

11,2 

1,0 

1.0 

— 

11,3 

1.0 

0 

— 

11,4 

1.0 

(9,8) 

500 

H,5 

1.0 

(5.8) 

1000 

11,6 

1.0 

(*.*) 

1500 

Resonance  spin  for  the  T108E40  is  l.l2°/ft.;  for  the  T316,  2.34°/ft. 

In  the  cases  1,3  through  1,6  and  11,4  through  11,6  it  is  seen  that 
from  any  two  conditions  the  third  may  he  derived.  The  derived  quantities 
are  indicated  by  parentheses. 
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Throughout  the  report  the  above  twelve  cases  will  be  referred  to  as 
1,1  .....  11,6,  wherein  the  Roman  numerals  designate  the  shell  and  the 
Arabic  numerals  the  particular  spin  history  under  consideration. 

The  pertinent  characteristics  of  the  two  shell  are  given  in  Table  2 
and  were  obtained  experimentally  from  previous  range  firings  or  wind 
tunnel  tests.1  Note  significant  differences  in  drag  coefficient,  damping 
rate,  pitching  frequency,  and  zero  spin  trim  angle,  X  ,  between  the  two 
shell.  For  example,  the  T108  trims  at  5-07  degrees  per  degree  of  tail 
asymmetry  while  the  T3l6  trims  at  1,83,  with  the  result  that  the  latter 
shell  is  much  less  sensitive  to  fin  deformation. 

The  following  section  describes  the  equations  of  motion  which  were 
programmed  for  the  analog  computer. 


An  aerodynamically  asymmetrical  shell  will  assume  an  equilibrium,  or 
steady  state,  condition  after  all  transients  have  died  out.  The  result¬ 
ing  angle  of  inclination  of  the  geometric  axis  to  the  trajectory  is 
the  trim  angle.  It  is,  in  general,  a  function  of  spin.  The  symbol 
1^,  used  throughout  the  report,  represents  the  value  of  the  trim  angle 

defined  for  zero  spin. 
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TABLE  2 


PERTINENT  CHARACTERISTICS  OF  90MM  SHELL 


T108E40 

T316 

L  (overall  length  in  cal.) 

10.00 

6.96 

d  ( diameter  in  f t . ) 

.2945 

.2945 

m  (weight  in  lbs.) 

15.60 

14.50 

c.  g.  (cal.  from  base) 

3.6 

2.98 

B  (lb. -in.2) 

727 

372.6 

A  (lb. -in.2) 

24.5 

22.74 

pd^/m 

1.25(10'^) 

1.32(10'^) 

cu 

1 

OJ 

.268 

.485 

5.34 

2.92 

w 

5.07 

1.83 

Period  of  yaw  ( f t . ) 

322.0 

154.0 

Pitching  frequency  (rad. /cal.) 

.00576 

.01203 

Damping  rate  (l/cal.)  =  H/2 

7.36(10'4) 

1.77(10-4) 

uo(ft./sec. ) 

2800 

2800 

C,  (l/cal.) 

3.24(lo"5) 

3.24(10"5) 

vq  (rad. /cal. ) 

9.76(10'5) 

5.l4(lO~5) 

kd 

O 

.14 

+  * 
.22-. 02 

V2 

7.07(10'4) 

1.98(10“^) 

If 

Kp  is  a  linear  function  of  Mach  number  over  the  Mach  number  range 
o 

considered  in  this  problem  (see  Figure  3).  The  small  Mach  number 
variation  of  the  T108E40  was  neglected. 
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TABLE  2  (continued) 

PERTINENT  CHARACTERISTICS  OF  90MM  SHELL 

T108E40 

T3l6 

1.19 

l.4l 

**  <=-8V 

•  95 

1.13 

€ 

kl 

1.05 

1.10 

Si 

-1.00 

-2.25 

Si 

-5.07 

-4.12 

€ 


EQUATIONS  OF  MOTION 

If  the  assumptions  of  zero  initial  yaw  and  zero  initial  change  in 

yaw  are  made,  and  if,  further,  the  very  small  Magnus  and  yaw  of  repose 

terms  are  neglected,  the  linearized  differential  equation  of  yawing  motion 

for  missiles  ha-ving  slight  configurational  asymmetry  may  be  written  in 

2 

the  following  form 


X "  +  (H  -  iv)  X'  -  MX  =  J  el9  ,  » 

e  (1) 

where  primes  denote  differentiation  with  respect  to  p,  the  axial  arc 
length  along  the  trajectory  in  calibers.  The  nomenclature  used  in  the 
above  equation  and  in  the  equations  to  follow  is  defined  in  the  Table  of 
Symbols  preceeding  the  Introduction. 


If  the  small  imaginary  term  in  the  expression  for  J 
and  if  the  substitution 


€ 


is  neglected, 


X 


t 


M 


X 

e 


(2) 
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is  made,  equation  (l)  may  be  simplified  by  dividing  through  by  the  trim 
angle,  A^.  Separating  the  new  variable  A/A^  =  £  into  its  horizontal  and 
vertical  components,  the  equations  assume  the  form  in  which  they  were 
programmed  for  the  analog  computer: 


£jj"  +  +  V  £v»  -  M(£h  -  cos  0)  =  0 


(3) 


V'  +  HV  '  V  '  M^v  "  sin  =  0 


00 


The  equation  of  rolling  motion  for  a  missile  with  rotational  symmetry 
may  be  expressed  in  terms  of  the  quantities  C  (damping)  and  S  (steady 
state  rolling  velocity)  obtained  by  the  present  roll  reduction  technique  . 


v'  =  C±  (S  -  v) 


(5) 


If  all  the  aerodynamic  forces  other  than  the  normal  force  are 

2  r 

neglected,  the  differential  equation  of  swerving  motion  becomes  1 ^ 


£ 

d 


x  +  iy" 


„  10  . 

M  e  a. 


=  JL  A  + 


igd 


K 


-2 


u 


(6) 


The  third  term  of  the  above  equation,  when  integrated  twice,  will  be 
designated  by  y„;  i.e. 


2  -2 

g  d  u  dp  dp 


(7) 


The  amount  of  swerve  caused  by  the  aerodynamic  forces  acting  on  the 
missile  is  small  in  comparison  with  the  drop  due  to  gravity,  y„.  Therefore, 
in  order  to  avoid  a  scaling  problem  on  the  analog  computer  which  would 
arise  if  equation  (6)  were  programmed  as  it  stands,  it  was  found  convenient 
to  separate  the  two  portions  of  the  equation  and  plot  the  resulting  motions. 
One  further  simplification  was  made  in  order  to  observe  the  effect  of  yaw 
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on  the  gravity  term,  viz.,  instead  of  y_,  a  quantity  Ay  was  computed 

U-  u 

such  that 


=  Sd 


-2  -2 

(u  -  u  )  dp  dp 


(8) 


where 


(9) 


Hi  -  j 

U  D0  (10) 

Ay^,  then,  represents  the  difference  between  the  gravity  drop  which 
occurs  when  the  yaw  drag  effect  is  considered  and  the  gravity  drop  which 
would  obtain  at  zero  yaw. 


Jn  is,  of  course,  a  function  of  Mach  number.  The  drag  function 
0 

determined  from  range  firings  of  the  T3l6  shell  is  given  in  Figure  3. 

The  effect  of  variation  in  J  with  Mach  number  for  the  Mach  number  range 

covered  by  the  T108E40  shell  in  flight  was  found  to  be  negligible  and 

hence  a  constant  value  of  Jn  was  employed  in  equations  (9)  and  (10) 

%' 

when  computing  the  motion  of  the  T108. 


Resolving  equation  (6)  (minus  the  gravity  term)  into  its  real  and 
imaginary  components,  the  following  quantities  are  defined: 


x 

dX 


T  M  cos  9 

,  JL  ?H  +  7  -2  j  | 

X  K2  l*t 


y  +  yc 

dT~ 


JL«V  + 


M  sin  9 

■  -2. : 


(11) 

(12) 


Equations  (3),  (4),  (5),  (8),  (9),  (10),  (ll),  and  (12)  were 
rammed  for  the  comj 
were  set  equal  to  zero. 


programmed  for  the  computer.  All  initial  conditions,  except  v  1  and  u  1 , 

o  o 
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table  3 


Case 

Spin  His  tory 

(Cf.  TABLE  2) 

Yaw  History 

\nax.#  0ccurs  at 

(For  X^».2°)  Z  (yards) 

■\[  x2  +  y^  (mils  at  2000  yds . ) 

for  X  *.2° 

€ 

v-*0 

AyG  (rails  at  2000  yds.) 

X£*.lt°  X£=.6°  X  -.8° 

X£*I/0° 

i,i 

Zero  constant 

spin 

1.0 

—  — 

38.26 

— 

— 

. 

— 

1,2 

Constant  spin 

of  1.9°/ft. 

.6 

— 

.19 

— — 

. 

— 

1,3 

Zero  steady 

•state  spin 

3.8 

1750 

.19 

-.016 

.15 

.38  .72 

1.18 

I,L 

Resonance  at 

$00  yards 

3.6 

65o 

»  1 

.083 

.32 

.72  1.30 

2.08 

1,5 

Resonance  at 

1000  yards 

3.7 

1200 

.19 

.025 

.25 

.63  1.22 

2.08 

1,6 

Resonance  at 

1500  yards 

3.8 

1750 

.19 

.016 

.15 

.38  .72 

1.18 

11,1 

Zero  constant 

spin 

.it 

— — 

5.66 

— - 

-  - 

—  ~ 

11,2 

Constant  spin 

of  1.0°/ft. 

.it 

— 

.05 

— - 

— 

— 

— 

H,3 

Zero  steady 

state  spin 

.it 

.05 

— 

— - 

. 

—  — 

II,  It 

Resonance  at 

500  yards 

It. 6 

700 

.05 

.68 

2.75 

7.87  - 

n,5 

Resonance  at 

1000  yards 

6.3 

1250 

.05 

.53 

3.50 

10.17  - 

11,6 

Resonance  at 

1500  yards 

7.7 

1850 

.05 

.22 

.82 

2.12  - 

— — 

fror  the  cases  in  which  resonance  does  not  occur,  the  values  given  represent  steady  state  yaw,  or  magnitude  of  yaw  after  transients 

\ 

have  damped  out.  Note  that  the  yaw  levels  attained  by  the  GELA  in  the  constant  spin  cases  are  verified  by  the  plots  of  Figure  10, 
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RESULTS 


The  quantitative  results  of  the  program  are  presented  graphically 
in  Figures  4  through  14  and  summarized  in  Table  3.  Some  of  the  qualita¬ 
tive  aspects  of  the  program  are  discussed  below. 


The  graphs  representing  the  various  spin  histories  (Figures  4  and  5) 
are  self-explanatory.  Note  that  the  rolling  velocity  of  the  T108E40 
must  decrease  in  order  to  reach  resonance  whereas  that  of  the  T3l6  must 
increase.  Thus  it  will  not  be  possible  to  directly  compare  the  cases 
wherein  resonance  occurs  at  500  yards  (1,4  and  11,4),  for  example,  since 
the  spin  rate  of  the  T108  passes  through  zero  (an  undesirable  feature, 
as  will  presently  be  shown)  while  that  of  the  T3l6  attains  a  steady  state 
value  of  9.8  degrees  per  foot. 


Graphs  of  the  horizontal  component  of  yaw  for  various  spin  histories 
are  given  in  Figures  6,  7,  8  and  9  for  a  particular  degree  of  "effective 
asymmetry*' .  The  effective  angle  of  asymmetry  of  a  shell  is  directly 
related  to  the  trim  angle  and  may  be  considered  to  be  the  asymmetry  which 
produces  a  trim  yaw.  This  asymmetry  may  be  present  due  to  a  variety  of 
causes,  such  as  cocked  fins,  squashed  fins,  deformed  boom,  etc.  and  is  not 
always  measurable  prior  to  firing.  It  can  be  inferred  from  experimental 
data  obtained  from  a  firing  by  observing  the  yawing  motion  at  trim  and 
using  the  following  relationship : 


MXt 


(13) 


A  value  of  X  =  .2  was  chosen  as  a  representative  figure  for  current 
manufacturing  tolerances  and  used  throughout  the  program.  However,  since 
appears  as  a  factor  in  the  equations  of  yawing  motion  [~cf.  equations 
(3)  and  (4)  and  recalling  that  i/l,  =  £  J  ,  the  magnitude  of  yaw  for  other 
degrees  of  asymmetry  may  readily  be  inferred  by  applying  the  appropriate 
scale  to  the  given  graphs. 


* 

Because  of  the  choice  of  very  special  initial  conditions  (\  =  \'  =0=0), 

the  vertical  component  of  yaw  is  merely  90°  out  of  phase  with  the 
horizontal  component  after  the  transients  have  died  out  and  hence  only 
one  component  was  recorded. 
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The  behavior  of  the  yawing  motion  undergoing  the  effects  of  resonance 
is  of  considerable  interest.  It  is  seen,  by  comparison  with  Figures  4  and 
5,  that  maximum  yaw  is  not  attained  at  the  same  time  that  the  rolling 
velocity  and  pitching  frequency  are  coincidental,  but  rather  at  some 
distance  approximately  20$  further  along  the  trajectory.  Moreover,  the 
maximum  amplitude  of  the  resonating  yawing  motion  is  a  function  of  the 
speed  with  which  the  spin  passes  through  the  resonance  region,  i.e.  the 
acceleration  of  the  roll  angle,  9.  For  the  cases  presented  in  this 
report  the  spin  velocity  is  greatest  where  resonant  spin  is  reached  at 
500  yards  and  least  for  resonance  at  1500  yards.  Furthermore,  both  the 
T108  and  the  T3l6  shell  attain  the  highest  yaws  when  resonant  spin  is 
reached  at  1500  yards.  Thus  the  more  slowly  spin  passes  through  the 
resonance  region  the  greater  is  the  maximum  amplitude  attained  by  the 
yawing  motion.  The  phenomenon,  moreover,  is  more  pronounced  in  the  case 
of  the  T316  shell;  this  aspect  will  be  explained  by  the  following 
discussion. 


For  constant  spin  a  steady  state  solution  of  the  form  K-.  e^v-^  can  be 

3 

assumed  for  equation  (l),  and  the  following  relation  obtained  for  the 
magnitude  of  K, : 


K, 


Mi  -  a/b)jh  -  k2  -2  j 

€  £  J 


(14) 


(l  -  A/B)v  -  ivH  +  M 


Neglecting  the  very  small  contribution  of  the  imaginary  term  in  the 

numerator  of  the  above  expression,  the  resonance  curves  for  the  two  shell 

o  -X- 

were  calculated  for  =  .2  and  are  shown  in  Figure  10.  The  abscissa 
of  each  point  on  the  curves  represents  the  level  of  yaw  that  would  be 
attained  for  the  particular  value  of  (constant)  spin  indicated  as  the 
ordinate  of  that  point.  If  these  two  curves  are  properly  used,  they  are 
of  considerable  value  in  interpreting  the  yawing  motions  obtained  by  the 
GEDA. 


In  mechanical  systems  in  general,  this  type  of  plot,  which  yields  the 
amplitude,  or  response,  of  the  system  to  an  external  force  of  any  given 
frequency,  is  known  as  a  response  curve. 
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It  is  seen  from  the  graphs  that  the  maximum  values  of  yaw  that  would 
be  attained  by  the  T108  and  the  T316  for  constant  resonant  spin  are 
approximately  4°  and  12°  respectively.  The  resonant  maximum  yaw  of  the 
T108  shell  is,  therefore,  approximately  (4/. 2)  x  X  ,  or  in  resonance, 
the  shell  suffers  amplification  of  its  inherent  effective  asymmetry 
twenty  times.  Resonant  amplification  of  the  asymmetry  of  the  TJl6 
shell  is  even  larger,  being  (12/  .2)  x  or  approximately  sixty.  The 
pronounced  difference  between  the  peaks  of  the  two  curves  is  due  to 
the  difference  in  damping  rates  of  the  two  shell,  that  of  the  T316  being 
only  about  one -fourth  the  damping  rate  of  the  T108  (see  Table  2). 


A  measure  of  a  shell's  sensitivity  to  its  trim  angle  under  resonance 

conditions  is  given  by  its  magnification  factor  (m.f.),  defined  as  the 

ratio  of  the  value  of  K*  I  at  resonance  spin  to  the  value  of  I  K, I  for 

'  >  1  131 

zero  spin  (i.e.lt).  Values  of  the  magnification  factor  for  the  T108  and 
the  T316  shell,  indicated  in  Figures  10  and  11,  are  approximately  4  and 
40  respectively. 


Figure  11  illustrates  the  behavior  of  resonating  yawing  motion  as  a 
function  of  roll  angle  acceleration.  It  consists  of  two  plots  ,  repre¬ 
senting  the  T108E40  and  the  T3l6  shell,  whose  coordinates  are  dimensionless 
and  defined  in  the  following  manner: 


Abscissa:  Maximum  Yaw  (v  variable J 


ax. 


K, 


at  Resonance 


K, 


Reson. 


Ordinate : 


A  z 


Period  of  Yaw 


=  Periods  of  Yaw  Near  Resonance 


where 


Az  -  z(v2)  -  z(vxs 


v 

v 


1 

2 


=  .95  v 

reson. 
=  1.05  v 

reson. 


(15) 


¥ . 

These  plots  were  suggested  by  G.  E.  Murphy. 
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In  other  words,  a  "region  of  resonance",  arbitrarily  chosen  to  lie  within 
-  5 1o  of  critical  resonance  spin,  was  marked  off  and  the  distance  along  the 
trajectory,  corresponding  to  that  region,  was  computed  for  each  of  the  six 
resonance  cases.  This  distance,  Az,  divided  by  the  period  of  yaw  (see 
Table  2)  becomes  the  parameter  against  which  the  percent  of  maximum  yaw 
for  constant  spin,  attained  along  the  recorded  trajectory,  was  plotted. 

(Note  that  maximum  yaw  does  not  necessarily  occur  within  the  "x-egion  of 
resonance"as  defined  in  the  present  discussion.)  The  results  of  the  two 
plots  are  smooth  curves  which  asymptotically  approach  unity.  Moreover, 
they  are  suggestive  of  a  family  of  curves  and,  since  it  seemed  of  some 
(at  least  academic)  interest,  a  side  program  was  set  up  to  test  that  theory. 
On  the  order  of  fifteen  additional  GEDA  runs  were  made  in  which  the  damping 
coefficient,  H,  in  the  yaw  equation  was  allowed  to  take  on  different  values. 
Results  showed  that  the  two  curves  given  in  Figure  11  are  indeed  members 
of  a  family  of  curves  having  the  magnification  factor  as  the  parameter. 

Turning,  now,  from  our  consideration  of  the  behavior  of  the  yawing 
motion  to  the  lateral  displacement  of  the  shell,  and  referring  once  again 
to  the  equations  which  were  programmed  for  the  analog  computer,  it  is  seen 
that  the  vector  sum  of  the  motions  (x,  y  and  A y_>)  pertaining  to  the  lateral 
displacement  obtained  by  the  simulation  process  is,  in  effect,  a  departure 
from  a  particle  trajectory  rather  than  the  actual  trajectory  itself,  x 
and  y  are  the  two  components  cf  the  swerving  motion  of  the  shell  due  to 
the  aerodynamic  forces  acting  upon  itj  Ay^  is  the  departure  of  the  shell 
from  a  particle  trajectory  due  to  increased  drag  alone.  In  virtue  of  very 
special  initial  conditions  (cf.  paragraph  following  equation  12),  x,  the 
horizontal  component,  is  so  small  for  most  of  the  runs  (y  is  greater  by  a 
factor  of  between  20  and  100 )  as  to  be  negligible  and  hence  the  graphical 
representation  of  this  motion  has  been  omitted  from  the  report.  In  actual 
practice,  the  orientation  of  the  departure  of  the  shell  could  occur  in  any 
plane  depending  upon  the  launching  conditions . 
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In  the  preceding  paragraph  it  was  stated  the  for  most  of  the  runs 
x  is  very  small.  Three  exceptions  occur:  Cases  1,1  and  11,1  in  which 


It 


the  spin  rate  is  identically  zero,  making  the  first  term  in  the  x  equation 

constant  after  the  transients  have  died  out  (the  second  term  is  constant 

hy  virtue  of  the  zero  spin  condition)  and  thereby  causing  the  solution  to 

be  almost  purely  parabolic;  and  Case  1,4,  in  which  the  spin  rate  passes 

rather  slowly  through  zero  at  approximately  4000  feet  from  the  beginning 

of  the  range.  Since,  for  constant  spin,  the  deflection  for  missiles  having 

2 

slight  configurational  asymmetry  is  proportional  to  l/v  ,  it  is  evident 
that  the  condition  of  zero  rolling  velocity  is  to  be  avoided. 

Figure  12  contains  plots  of  y  for  two  representative  cases.  It  was 
found  that  the  passage  of  the  missile  through  resonance  at  different  places 
along  the  trajectory  had  no  appreciable  effect  upon  x  or  y;  i.e.,  the 
radius  of  the  helix  caused  by  the  yawing  motion  is  negligible  in  comparison 
with  the  jump  due  to  asymmetry.  For  this  reason,  the  other  ten  cases  are 
not  recorded  since,  for  most  of  them,  the  amount  of  departure  is  the  same 
as  for  the  two  cases  shown.  Actually  y  for  the  zero  spin  cases  (1,1  and  II, l) 
is  identically  zero  in  virtue  of  the  specific  initial  conditions  of  the 
problem.  Case  1,4  is  also  exceptional  since  y  becomes  very  large  when  the 
spin  rate  passes  through  zero. 

All  of  this  information  is  tabulated  in  Table  3  under  the  column 


The  data  in  this  column  are  the  absolute  values  of 


the  sum  of  x  and  y  expressed  in  mils.  Note  that  the  values  for  the  T108E40 

are  on  the  order  of  25$>  higher  than  those  for  the  TJl6  because  of  the  fact 

that  the  T108  has  a  much  higher  degree  of  trim  yaw.  The  effect  of  the 

difference  in  trim  angle  between  the  two  shell  is  especially  pronounced 

in  the  zero  spin  case  wherein  it  has  been  shown  that  the  x  motion  is 

almost  purely  parabolic.  For  that  case  the  departure  of  the  T108  is  about 

seven  times  that  of  the  T3l6  at  a  target  of  2000  yards.  As  in  the  case 

of  the  yawing  motion,  the  values  given  are  for  an  effective  asymmetry  of 
o 

.2  ;  for  other  degrees  of  asymmetry  it  is  necessary  to  multiply  by  the 
appropriate  scale  factor. 
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The  remaining  data  obtained  in  the  program  are  concerned  with  the 
effect  due  to  yaw  drag  on  the  gravity  term  in  the  swerve  equation. 

Because  of  the  non-linear  nature  of  the  equation  for  y^,  simple  scaling 
of  yaw  effects  is  impossible,  and  computations  were  made  for  angles  of 
asymmetry  between  .2°  and  1°  rather  than  for  one  value  alone  as  in  the 
case  of  the  yaw  and  swerve  equations.  Since  fin  distortions  of  1°  and 
even  higher  have  actually  been  observed  on  firings  of  the  T108  production 
shell,  the  range  chosen  for  this  study  is  realistic.  It  is  felt,  however, 
that  by  its  very  nature  the  T3l6  is  not  likely  to  suffer  fin  damage  of 
very  large  magnitude.  Indeed,  because  of  the  extremely  high  level  of  yaw 
that  would  be  encountered,  assumption  of  as  much  as  1°  of  tail  malalignment 
of  the  T316  does  not  permit  the  use  of  the  Kelley-McShane  linearized 
equations  of  yawing  and  swerving  motion  which  are  based  on  the  assumption 
of  small  yaw.  For  these  reasons,  therefore,  the  study  of  the  yaw  drag 
effect  on  the  T316  was  limited  to  yaws  based  on  .6°  of  asymmetry  or  less. 

A  few  plots  of  the  AyQ  motion  (Figures  13  and  Ik)  are  presented  to  show 
the  type  of  non-linear  effect  due  to  yaw  drag  caused  by  an  increase  in  yaw 
level.  Note  that  the  amount  of  deflection  at  a  1000  yard  target  is 
practically  insignificant  in  all  cases.  However,  it  is  apparent  that  the 
deflection  at  2000  yards  is  definitely  a  function  of  the  place  of  occurrence 
of  resonance  along  the  trajectory  (See  also  Table  3).  If  early  resonance 
is  avoided,  and  fin  damage  remains  small,  the  yaw  drag  effect,  even  for 
the  T3l6  shell,  is  not  too  detrimental. 

CONCLUSION 

In  order  to  summarize  the  results  of  the  program  and  draw  certain 
conclusions  therefrom,  it  is  worthwhile  to  review  its  fundamental  objectives. 
In  the  introduction  an  attempt  was  made  to  review,  at  least  in  part,  the 
history  behind  the  T108E40  and  the  T3l6  shell.  With  this  background  in 
mind,  then,  the  problems  which  emerge  from  the  picture  are  mainly  three 
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in  number: 


(1)  Does  the  resonance  phenomenon  account  for  the  poor  firing 
accuracy  of  the  T108  shell? 

(2)  Is  the  ''time  of  resonance"  (place  wherein  resonance  occurs 
during  flight)  significant? 

(3)  Could  the  reason  for  the  good  accuracy  of  the  T316  shell 
be  the  fact  that  it  is  inherently  a  better  shell,  from 
aerodynamic  considerations,  than  the  T108? 

Resonance  undoubtedly  accounts  f  oi  some  of  the  poor  firing  accuracy 
observed  in  range  firings  but  its  significance  appears  to  be  less  than 
was  suspected  at  the  outset  of  this  piogram.  Dispersions  of  several  mils 
at  1000  yards,  determined  experimentally,  have  not  been  obtained  by  the 
analog  computer  study.  It  appears  that  the  observed  behavior  of  production 
lots  of  the  T108  shell  has  not  been  simulated  since,  from  this  study,  the 
shell  appears  to  fly  very  well  to  a  1000  yard  target.  At  the  time  of  this 
writing,  an  investigation  is  being  conducted  '  to  determine  the  effect  on 
trajectory  of  initial  yaw  and  yawing  velocity  in  an  effort  to  explain  the 
poor  firing  accuracy  of  the  T108  shell.  The  results  appear  promising. 

Question  two  is  quite  easily  resolved  from  a  study  of  the  results 
obtained  from  the  analog  computer.  It  is  seen  (Table  3)  that  the  amount 
of  displacement  due  to  the  aerodynamic  forces  acting  on  the  shell  is  not 
a  function  of  the  place  wherein  resonance  occurs  during  flight.  The 
maximum  amplitude  of  the  yawing  motion  does  depend  upon  the  "time  of 
resonance",  but  this  relation  is  important  only  when  considering  its 
effect  upon  the  total  trajectory.  That  is  to  say,  the  variation  of 
maximum  yaw  level  with  time  of  resonance  is  of  little  more  than  academic 
interest  as  far  as  the  yawing  motion  alone  is  concerned,  whereas  it  has 
considerable  significance  with  respect  to  the  amount  of  displacement  due 
to  yaw  drag.  Since  the  drag  coefficient  is  a  function  of  yaw  squared, 
the  relationship  of  Ay^  to  the  size  of  the  yawing  motion  is  non-linear 
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and  at  2000  yards  its  contribution  to  the  total  lateral  displacement  is 
of  considerable  magnitude.  Nevertheless,  at  1000  yards  the  departure  is 
almost  insignificant  and  hence  yaw  drag  cannot  be  the  explanation  for  the 
poor  firing  accuracy  of  the  production  model  T108  shell  whose  dispersions 
at  a  1000  yard  target  were  observed  to  be  on  the  order  of  one  mil. 

Unfortunately,  question  three  cannot  be  answered  unequivocally.  A 
comparison  of  the  aerodynamic  properties  of  the  two  shell  would  provide  an 
adequate  explanation  for  the  better  firing  accuracy  of  the  TJl6  if  we  had 
been  able  to  prove  that  resonance  were  really  the  cause  of  the  bad  behavior 
of  the  T108.  In  resonance,  the  T108  shell  suffers  amplification  of  its 
inherent  asymmetry  twenty  times.  The  tail  of  this  shell  is  relatively 
weak,  especially  if  it  is  made  of  aluminum  63ST;  it  is  known  to  suffer 
considerable  deformation  during  launching  on  at  least  50%  of  tne  rounds. 

The  resulting  effective  asymmetry  is  often  as  large  as  one  degree  or  even 
larger.  Moreover,  with  relatively  low  initial  spin  this  shell  can  get  into 
resonance  quite  early  along  its  trajectory. 

On  the  other  hand,  resonant  amplification  of  the  T3l6  shell  is  consider¬ 
ably  larger :  it  is  about  sixty  times  its  inherent  asymmetry,  however,  its 
resonant  frequency  is  approximately  2.J4  degrees  per  foot  and  starting  with 
zero  or  very  low  spin,  therefore,  it  requires  a  relatively  large  deformation 
of  the  tail  to  reach  resonant  spin  level  within  short  ranges.  Because 
this  shell  is  more  compact,  it  is  likely  to  launch  better,  i.e.  with  less 
deformation.  Hence  its  chances  of  reaching  resonant  spin  are  much  smaller 
than  those  of  the  T108.  A  comparison  of  the  resonance  curves  shown  in 
Figure  10  reveals  that  there  is  over  twice  as  large  a  region  of  spin  wherein 
the  T3l6  will  not  be  near  the  critical  resonance  value.  Therefore,  in  this 
sense  it  is  a  better  shell,  and  less  likely  to  experience  detrimental  effects 
due  to  resonance. 


The  results  of  this  program,  although  not  providing  the  complete 
solution  to  the  dispersion  problem  of  the  T108  shell,  nevertheless  shed 
sufficient  light  on  the  situation  so  as  to  indicate  worthwhile  paths  for 
future  investigations.  ,  , 

J.  M.  SCHMIDT 
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FIGURES 


1.  Photograph  of  T108E40 

2 .  Photograph  of  T3l6 

3.  Drag  Function  of  T3l6  (K^  vs.  Mach  Number) 
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8.  vs  Z  for  T316  (cases  11,1;  II, ,2;  11,3;  and  11,5) 
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10.  Resonance  Curves  (K^  vs.  Spin) 
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APPENDIX:  COMPUTER  CHECK  SOLUTIONS 

Several  special  runs  were  made  on  the  computer  for  the  purpose  of 
checking  the  validity  of  the  solutions  being  obtained.  They  were  arranged 
so  as  to  check,  as  nearly  as  possible,  all  of  the  parts  of  the  circuit. 


When  the  rolling  velocity  of  the  missile  is  constant,  the  solution 
to  the  differential  equations  of  yawing  and  swerving  motion  may  be 
witten  exactly  in  terms  of  simple  functions,  thereby  permitting  a  fairly 
simple  check  of  the  computer  solution  for  constant  spin.  These  equations 
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the  gravity  term  having  been  omitted  from  the  present  discussion.  With 
the  assumption  of  zero  initial  yaw  and  zero  initial  change  in  yaw. 
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The  spin,  v,  may  be  chosen  arbitrarily  for  this  computation,  while 
and  the  damping  rates  and  turning  rates,  respectively,  of  the 
two  epicyclic  arms,  are  specified  by  the  equations 


$1  .  $2  =  -  M 
$1  +^2  =  “  H  +  1  V 


(24) 

(25) 


X.^.,  the  trim  angle,  or  magnitude  of  yaw  due  to  asymmetry  after  the 
transients  have  died  out,  is  given  approximately  by 


J..  X 

M  € 


M 


(26) 


Neglecting  the  very  small  imaginary  term  in  equation  (20),  K  may  then 

3 

be  rewritten  in  terms  of  the  trim  angle,  X. : 

T? 


k,  „ 

3  _  _ _ 

\  '  (fx  -$3)  ($2  -§j) 


(27) 


If  the  constant  and  linear  terms  of  equation  (17)  are  designated 
as  P/dX^.  and  Q/dX^,  respectively,  where  P  and  Q,  are  complex,  then 


i=l 


K. 


+ 


M 


i  -2  2  | 

k2  V  k 


M 


h2  "2  v 


(28) 

(29) 
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Since  the  remaining  terms  of  the  swerve  equation  are  oscillatory  in 
nature,  it  is  necessary  to  compute  only  (P1  +  Q^p)  and  (Pg  +  Q^p)  in  order 
to  determine  the  horizontal  (x)  and  vertical  (y)  displacements  at  some 
given  distance,  p,  in  calibers  along  the  trajectory.  With  the  following 
data  as  input, 

v  =  5*1^-  x  10  *  rad. /cal. 

H  =  1.4725  x  10"5 

M  =  -.52964  x  10"^ 

\t  =  1.104 

d  =  .2945  ft. 

the  necessary  quantities  for  solving  equations  (28)  and  (29)  were 
computed,  with  the  result  that  at  p  =  20575.5  cal.  (2000  yards) 

HMD  COMPUTATION  COMPUTER  SOLUTION 

x  .028  ft.  .05  ft. 

y  2.175  ft.  2.15  ft. 


In  order  to  obtain  a  check  solution  of  the  gravity  term  in  the 

equation  of  swerving  motion,  a  constant  drag  coefficient  was  assumed, 

thus  permitting  the  direct  formal  integration  of  the  expression  for  y_ . 

G 

Thus,  for  constant  J^, 


yG  =  gd 


rP  op  -2  2  -2 

J  J  u  dp  dp  =  gd  UQ  « 


e  ^  -2Jdp  -  1 

<77 


(30) 


where  u  =  u  e 
o 
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With  the  following  data  as  input. 


u  =  2800  ft. /sec. 
o  ' 

g  =  32.16  ft. /sec. 2 
d  =  .2945  ft. 

JD  =  4.848  x  10"5 
P  =  20373.5  cal. 

a  hand  computation  of  the  quantity  y^  was  made  and  compared  with  the  analog 
computer  solution  with  the  result  that 

HAND  COMPUTATION  COMPUTER  SOLUTION 

u  1042.8  ft. /sec.  1040.3  ft. /sec. 

Yq  160.12  ft.  161.00  ft. 

A  similar  procedure  was  followed  to  check  the  inclusion  of  the  yaw 
drag  term.  A  constant  average  value  of  yaw  squared  of  4  square  degrees 
was  assumed,  and  the  velocity  equation  was  permitted  to  take  the  form 

-JD» 

u  =  u  e 
o 

2  -7 

where  2  ^  2  ~  ^*^08  x  10  .A  non-linear  unit  of 

o  6  6 

o 

the  analog  computer  was  used  to  compute  J  _  5  ,  even  though  for  this 

V 

computation  all  the  quantities  are  constant,  in  order  to  simulate  as 
closely  as  possible  the  actual  requirements  of  the  problem  when  the 
yawing  motion  is  not  constant.  Results  again  indicated  a  fairly  high 
degree  of  accuracy  for  the  analog  computer  solution. 

HAND  COMPUTATION  COMPUTER  SOLUTION 

yQ  156.9  ft. 


157.2  ft. 
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